function [CC] = contactCentral(colsOutToMid)
% colsOutToMid = matrix (Nsamples,4)
%   each column in matrix represents a column in the palate,
%   col 1 being the lateral-most, and col 4 being the center-most (only
%   looking at one side per use of the function)

% Can operate on entire matrix of values, where rows represent different
% samples

CC = log((1*(colsOutToMid(:,1)/14)...
    +17*(colsOutToMid(:,2)/16)...
    +289*(colsOutToMid(:,3)/16)...
    +4913*(colsOutToMid(:,4)/16))+1)/log(5220+1);

end